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possibly strong solar flare or a very strong cosmic ray event of unknown origin that occurred 



ABSTRACT 

Aims. Miyake et al. (henceforth M12) recently reported, based on '''C data, an extreme cosmic 
OO ', event in about AD775. Using a simple model, M12 claimed that the event was too strong to be 

' caused by a solar flare within the standard theory. This implied a new paradigm of either an im- 

(N ■ 
O 

CO ' around AD775. However, as we show, the strength of the event was significantly overestimated 

by M12. Several subsequent works have attempted to find a possible exotic source for such an 

event, including a giant cometary impact upon the Sun or a gamma-ray burst, but they are all 

, based on incorrect estimates by Ml 2. We revisit this event with analysis of new datasets and con- 

$H < 

, sistent theoretical modelling. 

Methods. We verified the experimental result for the AD775 cosmic ray event using indepen- 
dent datasets including '"Be series and newly measured '''C annual data. We surveyed available 
historical chronicles for astronomical observations for the period around the AD770s to identify 
potential sightings of aurorae borealis and supemovae. We interpreted the '"'C measurements us- 
ing an appropriate carbon cycle model. 

Results. We show that: (1) The reality of the AD775 event is confirmed by new measurements 
of '''C in German oak; (2) by using an inappropriate carbon cycle model, M12 strongly overes- 
timated the event's strength; (3) The revised magnitude of the event (the global ''*C production 
2=(1.1-1.5)-10* atoms/cm^) is consistent with different independent datasets ('''C, '"Be, ^^*C1) 
and can be associated with a strong, but not inexplicably strong, solar energetic particle event (or 
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a sequence of events), and provides the first definite evidence for an event of this magnitude (the 
fluence > 30 MeV was about 4.5 ■ 10'" cm"') in multiple datasets; (4) This interpretation is in 
agreement with increased auroral activity identified in historical chronicles. 
Conclusions. The results point to the likely solar origin of the event, which is now identified as 
the greatest solar event on a multi-millennial time scale, placing a strong observational constraint 
on the theory of explosive energy releases on the Sun and cool stars. 

Key words. Sumactivity - Sun:flares 



1. Introduction 



It is of particular interest, both from astrophysical and societal points of view, to understand the full 
spectrum of severity of extreme events in the Earth's radiation environment. In particular, knowing 
the maximum possible energy of solar energetic particle (SE P) events and the frequency of their 
occurrence is of great importance for solar and stellar physics jHudsonl201ol) . The history of direct 
solar observations is relatively short, spanning only several d ecades and provide s insufficient statis- 
tics on the occurrence rate of the most energetic SEP events ( Smart et al.ll2006 ). Instead, long-term 
records of cosmogenic isotopes, like and '"Be stored in terrestrial, m eteoritic, or lunar archives, 
can re veal the history of such events over a greatly extended time-span ( Usoski JI2OO8 



20121) . While the statistics of extreme SEP events remains unclear to some extent ( Hudson 



progress in this field have been recently made (ISchriiver et al 



2012 : 



Beer et al 



2010) 



Usoskin & Kovaltsov 



2012h 



allowing the occurrence frequency of such events to be better estimated by joint analysis of differ- 
ent cosmogenic isotope data. 

A particularly exciting result has been recently published by 



Mivake et al 



(120121 henceforth 



Ml 2) who found a significant enhancement of about 1.5% (15 permill) of '"'C content measured 
in Japanese cedars around AD775. Using a basic four-box model of the carbon cycle, M12 esti- 
mated the corresponding absolute global production for the event as 6 • 10** atoms/cm^ or 19 
atoms/cm^/sec averaged over a year This is an order of magnitude greater than the average '"^C 
production rate due to galactic cosmic rays (GCR), which is estimated to be 1.6-2 atoms/cm^/sec 
( Kovaltsov et alj|2012 , and references therein) for the pre-industrial era. When translating the pro- 
duction rate into the flux of cosmic rays or the energy of the source, candidates being either a giant 
solar eruption or a nearby supernova, M12 concluded that it was much too high, implying a su dden 
strong cosmic ray event of unknown origin. In an attempt to resolve the situation lAUenl (120121) pro- 
posed "a supernova largely hidden behind a dust cloud . . . The resulting supernov a remnant would be 
invisib le." However as discussed below, this interpretation is unlikely. Recently 



Melott & Thomas 



(120 12I) recalculated the energy of a possible coronal mass ejection (CME) related to such a burst of 
SEPs and found a lower value. But it is still too energetic to correspond to a r ealistic solar eruptive 
event, leaving open the problem of the possible source of the AD775 event. 



Eichler & Mordecai 



(120 12h tried to link the event source to an impact of a cometar y body upon the Sun producing a 



Hambarvan & Neuhauserl (120131) 



shock in the corona to sufficiently accelerate energetic particles, 
proposed a gamma-ray burst as a possible source. Thus, a number of exotic sources have been 
proposed, because the most plausible one, the SEP event, looks unrealistically energetic given the 
strength of the event as estimated by Ml 2. Instead, the M1 2 result disagrees wi th data on another 



2008). The two cosmo- 



cosmogenic radionuclide, '"Be in the Dome Fuji ice core (iHoriuchi et al. 
genie isotopes ('"^C and '"Be) are formed as sub-products of the same process of a nucleonic- 
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Fig. 1. Time profiles of the measured A'^^C content in Japanese cedar (IMivake et al. 
German oak (this work) for the period around AD775. Smooth black and grey hnes depict a family 
of best fit A'^C profiles, calculated using a family of realistic carbon cycle models for an instanta- 
neous injection of '"^C into the stratosphere. 
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Fig. 2. A scheme of the carbon cycle models used here with aiTows depicting the carbon exchange 
between different reservoirs. Big arrows denote production Q of '"*C in the atmosphere. 



electromagnetic -muon cascade c aused by energetic cosmic rays or y-quanta in the Earth's at- 
mosphere (e.g.. iBeer et al.l 120121) . Thus, it is hardly possible to obtai n a tenfold increase in the 



20061). This 



annual '"^C production without a corresponding increase in "'Be data (lUsoskin et al. 
indicates an inconsistency i n the scenario proposed by M12. This event was also analyzed by 
Usoskin & Kovaltsovl ( 120121) who found the strength of the event much lower than did Ml 2. Here 
we revisit the AD775 event with analysis of new datasets and consistent theoretical modelling. 



2. ^"^C in German oak 

In order to verify the existence and strength of the claimed AD775 event, annual samples from 
a German oak (tree S teinbach 91 from the river Main), a part of the German Oak Chronology 
dFriedrich et al were measured independently in facilities of Mannheim, Germany (MAMS) 

and the ETH Zurich, Switzerland (ETH) for the period AD770-7 80. The data (Fig.[Tl) are presented 
here for the first time. One can see that they confirm the increase of production around AD775 
found by M12. Despite a possible slight vertical offset explainable by local peculiarities, the mag- 
nitude and timing of the increase are in full agreement. This implies that the event was global and 
was caused by the enhanced production of ^'^C. 



3. Analysis of ^^C data 

Only the relative content of '"*C is measured in the atmosphere. In order to evaluate the '"^C pro- 
duction, and hence the strength of the event, one needs to run a carbon cycle model and to compare 
the relative increase of A'^^C with the background concentration of '"*C in the troposphere, from 
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where it is absorbed by living trees. Here we use a family of th ree carbon cyc le models, includ- 
ing the main ca rbon reservoirs (Fig. |2): (1) a five-box model (pormanll2004 page 703); (2) a 



five-bo x model ( Dergachev & Vekslei 



model (lOeschger et al.lll975 



1991 



Damon & Peristvkh 



2004 . and (3) a box-diff^usion 



Siegenthaler et al.lll980l) . For better modeling of the SEP signal we 



also divided the atmosphere into stratosphere and troposphere, with a stratospheric carbo n capacity 



of 15% (cf. iMivake et alj|2012l) and a residence time of two years (IDamon et al 



1978). A model 



(IKovaltsov et alJl2012r) was used to simulate the C production by SEPs with a hard energy spec- 
trum, assuming 70% of ''^C produced by SEPs in the stratosphere and 30% in the troposphere, 
globally. These carbon cycle models diff'er slightly in exchange times between reservoirs and in 
dealing with the ocean, but they yield very similar results for the expected signal (Fig.[T]|. The 
background ^ '^C level was taken as the mean pre-industrial production rate of 1 .6 atoms/cm^/sec 
JoosiJbooiL The best fit by the weighted least-square method between the three datasets and the 
three models yields a net '"^C production of (1.3 + 0.2) ■ 10^ atoms/cm^. This is a factor of about 
5 smaller than the M12 model (6 ■ 10^ atoms/cm^). M12 used a different carbon cycle model, ne- 
glecting the deep ocean (Fig.|2]i, which is the greatest reservoir containing 92-95% of all carbon. 
Thus, M12 greatly overestimated the background '"^C concentration in the troposphere. This led 
to an eiTor when translating the relative increase of A''*C into the '^C production and so the event 
strength. We conclude that M12 overestimated the production by a factor of 4-6 because of the 
inappropriate model. Accordingly, all energy/particle flux estimates based on the numbers given 
by M12 are incorrect by the same factor 

We propose that the total '"*C production for the AD775 event coiTesponds to a SEP fluence 
(> 30 MeV) off 30 -4.5 - 10'" cm^^ with a hard spectrum as per the SEP event of 23 Feb. 1956 
(cf. 



Usoskin & Kovaltsov 



20121) . Since C is produced by more energetic particles, the value of 



FiQ depends on the assumed spectrum. For example, by assuming a very soft SEP spectrum as per 
the event of Aug. 1972, one would obtain the F30 fluence 1 .8 ■ 10" cm"^. More robust is the fluence 
of SEP with energy > 200 MeV, which is (5.5-1 1) ■ 10^ cm"^ for our scenario irrespective of the 
assumed SEP energy spectrum. 

In addition to the annual measurements in individual trees, discussed above, we use a five- 



year averaged INTCAL09 global '"*C series (IReimer et al.ll2009l) (Fig.|3]\). The proposed scenario 
is consistent with the data (x^/DoF - 0.68). A t ime shift of a few years can be explained by the 
filtering of the raw data in the INTCAL09 series ( Hogg et al.|[2009 ). 



4. "'Be in ice cores 



Using the above scenario, we calculate d the expected increase in '"Be by applying the production 
model of IKovaltsov & Usoskiii (201C) an d assuming an instan t injection and intermediate atmo- 
spheric mixing of '"Be ( McCracke JI2OO4 : Vonmoos et al.ll2006 ). Note that the ratio of the produc- 



tion of '"*C and '"Be by a SEP event is almost independent of the assumption of the SEP spectrum 



Dome Fuji series (IHoriuchi et al 



(lUsoskin et alJl200 6 |). Here we analvze t wo B e series measured in ice cores. One is the Antarctic 



2007 



2008h (Fig. |3p). The proposed scenario perfectly fits the 



data by magnitude, but the peak is delayed by several years and formally appears after AD780. 
However, since the dating of the '"Be series is related to the resolution of the stratigraphic intervals 
sampled and is dependent on an age-depth model between tie points, dating uncertainties of several 
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Fig. 3. Time profiles of the measured cosmogenic isotopes (histograms with eiTor bars) for the 
period around AD775. SoHd dots represent the enhancement expected in each data series assuming 
the AD775 event scenario discuss ed here (^'^C product ion g = 1 .3 ■ 10^ at/cm^). A) INTCAL09 5-yr 
samples glob al atmospheric A'^C (IReimer et al.l2009l) : B) Quasi-decadal "'Be content in the Dome 



Fuji ice core jHoriuchi et al.l2008l) . our model data is shifted by 5 year s to match the observed data; 



C) Quasi-decadal '"Be content in the GRIP ice core (lYiou et alJll997h . 



2000 



Horiuchi et al. 



20071). Thus, our proposed 



years are possible in the ice-core datasets (iBeerI t 
scenario is totally consistent with the ^°Be data measured in the Antarctic Dome F uji ice core. The 
other dataset is the Greenland GRIP series (lYiou et al.ll997 ; Vonmoos et al LoOfii (Fig.Ill::). There 



is a weak increase during the AD770s but it is not pronounced. The proposed scenario yields an 
expected peak that is higher (by about 2cr) than the observed peak. Thus, the GRIP series is not 
fully consistent with our scenario and the other data series, but the existence of the peak cannot be 
excluded at the 5% level. Instead, the corresponding ^^Cl peak in the GRIP corj^ agrees with the 
proposed scenario. 



5. Historical observations of tlie Aurora Borealis 



We surveyed publishe d aurora catalogues from oriental chronicles at low latitudes 
and 



Keimatsu 



1973) 



Yau et al 



(119951) cite credible observations from Shanxi Province, China, in AD770 (twice), 
AD773, and AD775. The next near est observatio ns are at AD767 and AD786. We also survey cat- 
alogues from occidental chronicles. iLinkI (1 1 962h cites the "red cross" in the sky dated AD773/774 
in different manuscripts of the Anglo-Saxon Chronicle (England), and "inflamed shields" in the 
sky (Germany, AD776). The next nearest occidental observations are at AD765 and AD786. At 
this time the Bible was a key reference in interpreting natural phenomena, explaining the cryptic 
reporting of aurorae (e.g., the above cited "red cross", also translated as "red sign of Christ" by 
Swanton 2000). In a new survey of occidental chronicles, we identified probable aurorae in AD772 
("fire from heaven", Ireland) and in an AD773 apparition interpreted by Christians as riders on 
white horses (Germany). Furt her new Irish observations are dated AD765 and AD786, and are 



found in the Annals of Ulster (iMac Airt & Mac Niocaill 



19831). and fl ie above-cited AD773 report 



from Germany in the Royal Frankish Annals (IScholz & Rogersll970b . We surveyed historical nova 



Estimate by JB based on unpubhshed preliminary data. 
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Allen 



catalog ues (e.g. JXu et al.ll2000h but found no credible report of supernovae in the AD770s. 
(120 12I) interprets the red cross (Anglo-Saxon Chronicle) as an exotic nearby supernova with an un- 
observable remnant, but we interpret this as an aurora. This is supporte d by a report tha t the same 
year "snakes... seen extraordinarily in the land of the S outh-Saxons" dSwantonluOOOf) . Serpents 
often feature in descriptions of aurorae jDair01molll98oi . reflecting the sinuous movement of au- 
roral sttuctures.^^fealso note the uncertain dating of these events in the Anglo Saxon Chronicle, 
with ISwantonI j2000l) re-dating them to AD776, i.e., after the onset of the AD775 event. We do 
not directly associate any particular aurora with the '*C event, but a distinct cluster of aurorae be- 
tween AD770 and AD776 suggests a high solar activity level around AD775. With the next nearest 
observations around AD765-767 and AD786, this suggests an 1 1-year cyclicity. 



6. Discussion and Conclusions 



The existence of the AD775 event is confirmed using a larger dataset, including two new annual 
'"'C series from German oak (Fig.[T}, but the event's interpretation by M12 is found to be incorrect. 
Because of an inappropriate carbon cycle model, M12 strongly overestimated (by a factor of 4-6) 
the strength of the event. Here we re-conducted the analysis, using more realistic models. The con- 
sequent event-integrated '"*C production rate is (1.1-1.5)-10^ atoms/cm^. We have verified that this 
value is in agreement with all the considered data series, including '"Be records in polar ice, yield- 
ing a consistent view corresponding to a strong SEP event (or a sequence of events) with a hard en 



ergx 



spectrum, 25-50 times stronger than the SPE event of 23 Feb. 1956 (cf. 



Usoskin & Kovaltsov 



2010) 



20121) . Such an event, while very strong, is not impossible for the solar dynamo (iHudson 
Moreover, it can be a sequence of SEP events as, for example, happened in the Autumn of 1989, 
thus further reducing the severity of individual events. This is corroborated by the steep tail of 
the SEP event fluence distribution dSchriiver et al. 2012 ). Several po tential candidates for similar 
events have previously been identified (lUsoskin & Kovaltsovll2012h . but observed in single data 
series only, thus leaving room for a possible terrestrial origin (e.g. regional climate excursion). The 
AD775 event is the only one consistently observed in several independent datasets, thus providing 
the first unequivocal observational evidence of such a strong SEP event on a multi-millennial time 
scale. This places a strong observational constraint on the upper limits of solar eruptive events, 
which is important for solar and more broadly stellar physics. 

In conclusion, by correcting the M12 model, by providing new independent '"^C data, and by 
surveying available historical chronicles and published aurora catalogues, we revisited the AD775 
event to demonstrate that it can likely be attributed to a strong solar SEP event. We show that: 

* The existence of the AD775 event is confirmed by new measurements of ''^C in German oak 
and by the existing '"Be data from polar ice cores. 



Mivake et al 



(120 12h overestimated the event's strength by a factor of 4-6. This di- 



rectly affects subsequent works based on this incorrect e stimate (e.g.. 



Melott & Thomas 



Eichler & Mordecaill2012 



Hambarvan & Neuhauseiil2013l) . 



2012 : 



* The revised event is consistent with different independent datasets and is associated with a 
strong, but not inexplicably strong SEP event (or sequence of smaller events), providing the first 
definite evidence for a SEP event of this magnitude from multiple datasets. 
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* This interpretation is in agreement with enhanced auroral sightings reported in historical 
chronicles for the period. 
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